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We review some formal aspects of cosmological solutions in closed string theory with 

duality symmetric "matter" emphasizing the necessity to account for the dilaton dynamics 

for a proper incorporation of duality. We consider two models : when the matter action is 

the classical action of the fields corresponding to momentum and winding modes and when 

the matter action is represented by the quantum vacuum energy of the string compactified 
• i— i . 

on a torus. Assuming that effective vacuum energy is positive one finds that in both cases 
the scale factor undergoes oscillations from maximal to minimal values with the amplitude 
of oscillations decreasing to zero or increasing to infinity depending on whether the effective 
coupling (dilaton field) decreases or increases with time. The contribution of the wind- 
ing modes to the classical action prevents infinite expansion. Duality is "spontaneously 
broken" on a solution with generic initial conditions. 
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1. Introduction 

Target space duality is known to be one of fundamental properties of string theory [1] . 
Being in some sense a symmetry between large and small scales it may play an important 
role in string cosmology [2] if the contributions of the winding modes are included [2-4]. 
In order for the string theory cosmological equations to be invariant under the "non- 
static" generalisation of duality (or " a model duality") [5] it is necessary to account for 
the dependence on the dilaton field which transforms under the duality [5-10]. It is the 
combined metric-dilaton system of equations that is invariant under the duality if the 
"matter" action is duality symmetric (i.e. contains contributions of both momentum and 
winding modes) [11]. 

In this paper we shall discuss some formal aspects of "cosmological" (i.e. time depen- 
dent) solutions in the case of duality symmetric "matter" with one of the aims being to 
illustrate the effect of the winding mode contributions on the evolution of the scale factor. 
We shall not consider "realistic" cosmological scenarios as this was already done in [12]. 
The present paper is a review and an extension of certain aspects discussed in [12]. 

Our choice of the matter actions will be based on the assumption of "adiabaticity" 
understood in the following sense. Let the space-time be a product of a time line and an N 
dimensional torus with radii a; = e Ai (which we shall usually take to be equal). We shall 
assume that the matter action can be represented by the classical action of the modes of 
the string theory compactified on the torus with Oj replaced by functions of time. The 
resulting action will be invariant under the duality transformation which inverts the radii, 
shifts the dilaton and interchanges momentum and winding modes. The corresponding 
equations will describe evolution of time dependent perturbations of the torus vacuum 
(from a different point of view these equations may be considered as a generalisation of 
the Kosterlitz-Thouless type renormalisation group equations [13]). 

We shall also consider the case when the matter action is represented by the (zero 
temperature) vacuum energy of the gas of momentum and winding modes (i.e. by the 
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partition function of the string compactified on the torus) which is automatically invariant 
under the duality transformation of a\ and the dilaton. Though the vacuum energy vanishes 
for the superstring a formal analysis of this case is useful since the resulting equations are 
similar to the system one finds in a more realistic finite temperature case [12]. 

We shall try not to specify the number N of space dimensions (i.e. keep the value 
of the central charge c ~ D — 26(= N — 25) arbitrary ) so that some of our results may 
apply, e.g., to the case of D = 2 string theory ( cosmo logical solutions with non-zero c in 
the matter-free case were discussed in [14-15]). 

We shall find that both in the "classical" and "quantum" matter cases the scale factor 
(the common radius of the torus) oscillates between finite maximal and minimal values . 
If the effective dilaton coupling is decreasing with time the amplitude of oscillations is also 
decreasing and the asymptotic value of the radius is the Planck one (\/c7). Though the 
duality is spontaneously broken by initial conditions it is thus restored asymptotically. The 
role of the winding mode contributions is to prevent an infinite expansion in the "classical" 
matter case [12] (and, modulo the issue of interpretation, to prevent a contraction to zero 
in the "quantum" matter case). 

We shall describe the expansion/contraction in terms of the original unrescaled ("a 
model") metric. Though the effective gravitational constant is then time dependent, this 
is the metric strings directly interact with and hence the one measured by stringy "rods" . 
Similar point of view seems to be expressed in refs. 16-17 where the question of which 
metric should be used to measure the expansion in string cosmology was discussed. The 
use of "a model " metric is particularly natural in the duality symmetric case since it is 
the unrescaled metric that has a simple transformation law under the duality. We shall 
assume that the dilaton does not have a non-perturbative potential, i.e. is "massless" . 

In Sect. 2 we shall present the basic "cosmo logical" equations for the two choices of 
the matter sources. In the case of the classical matter we shall consider a "mechanical" 
interpretation of the resulting system of equations (with the central charge playing the role 
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of an "energy" ) and also note a correspondence with the renormalisation group equations 
(for a related discussion of the RG equations see [18]). 

In Sects. 3 and 4 we shall study the behaviour of the solutions in the classical and 
quantum matter cases respectively. The qualitative analysis shows that the radius is oscil- 
lating with time between its extremal values with the amplitude of oscillations increasing 
to infinity or decreasing to zero depending on initial conditions for the dilaton. 

Sect. 5 contains some concluding remarks. We shall briefly discuss the "anisotropic" 
case in which there are "large" and "small" (internal) radii . 

In Appendices A and B we shall explicitly solve the basic "classical" and "quantum" 
systems of equations in the asymptotic region of large radius and demonstrate the existence 
of a maximal radius of expansion. 

2. Basic equations 

1. We shall consider only the leading order terms in the low energy expansion of the 
tree level gravitational effective action of the closed bosonic string theory [19] (we shall 
absorb the gravitational coupling constant into <j> and often set a'=l) 



The corresponding equations for the gravity plus matter action S = Sq + S m are 




(2.1) 



c = 



3>- 26 >' 




V 



(2.2) 



C = c + R + AD 2 4>-A(d(t)) 2 



! e 2 ^ 6S, 



m 



(2.3) 



2 V^G 6<f> 



We shall consider the following cosmological background 



N 




(2.4) 



a». 



= e *i(t) ; <j> = <j>{t) , N = D-l . 



For a proper interpretation of the matter action we shall discuss it is necessary to assume 
that at least some of coordinates periodic. 

It is useful to introduce the "shifted" (and rescaled by 2) dilaton field <p 

N 

V = 2(j)-^\ i , yf^G e~ 2(t> = e~ v , (2.5) 
i=i 

which is invariant under the duality transformation [6-10] 

Aj = -Xi , 4> = - \i , Cp = ip (2.6) 

(here i = 1, N ; more general duality transformations are obtained by combining these 
basic ones ). We are assuming that all the dimensions are compact. In general, the effective 
coupling is represented by the dilaton (2.5) where only the logarithm of the volume of the 
compact part of the space is subtracted out. Re- written in terms of Aj and <p the action 
and equations (2. 3), (2. 2) take the form (we drop a constant factor of integral over 

r , N 

S=- dt e^y^G^ [c-G 00 J2^ + G°V ] + S m [G 00 ,\ ,<p] , (2.7) 
J i=i 

N 

S = - dte-^[c + J2'^-^} , (2-8) 

N 

\-^\ = -\e^ , (2.10) 

^-E\ 2 = -e^-^). (2.1D 

Similar set of equations was discussed in [10-12]. Eqs.(2.9)-(2.11) follow from the action 
(2.7) (after the variation Goo is set equal to -1). If one stars from the action (2.8) (where 
Goo =~1) one should add eq.(2.9) as an additional "zero energy" constraint. If the original 



matter action is generally covariant the time derivative of (2.9) vanishes identically once 
(2.10) and (2.11) are satisfied. The examples of matter actions «S m [Goo> Aj, <p] we shall 
consider will be invariant under reparametrisations of time so that the identity relating 
the derivative of (2.9) to (2. 10), (2. 11) will be satisfied. Hence one of eqs.(2.10),(2.11) is 



integration constant which appears after integrating once eqs.(2.10),(2.11) is not arbitrary 
but is proportional to c. 

Eqs.(2.7)-(2.11) are invariant under the duality transformation (2.6) provided the 
matter action is duality symmetric. Given a solution of (2.9)-(2.11) we can then generate 
other solutions by applying duality transformations. 

Let us note that in the duality non-invariant case when the matter action is given 
by the dilaton potential V(4>) , i.e. (we shall assume the isotropic case when all radii are 
equal) 



2. As was noted in the Introduction, one of our aims is to understand the effect of the 
winding modes on cosmological evolution. Let us first consider the case when the matter 
action is represented by the classical action of the fields corresponding to momentum and 
winding modes of string theory compactified on a torus. For simplicity we shall consider 
just two complex scalar fields representing the momentum and winding modes of the 
tachyon field (xi are the "dual" coordinates ; a' = 1) 



redundant and can be dropped (if A ^ , ip ^ 0). The meaning of eq.(2.9) is that the 




NX 



the system of equations (2.9)-(2.11) reduces to 



c - NX 2 + (p 2 = e^V , 





(2.12) 



, N 

S m = \ dte^^G^[G 00 \i;\ 2 + (£ TO J e - 2A * -4)|V>| 2 
J i=i 

■ 2 N 

+ G 00 $\ + QTm 2 e 2A *-4)H 2 + 2^,V0] , (2.13) 
i=i 

where the potential U starts with quartic terms (cubic interactions involve other modes as 
well). To desribe the case when tp and ifi represent higher level scalar modes the tachyonic 
mass term -4 should be replaced by m 2 , = 4(n — 1) , n = 0,1,... . Variation of the 
total action (2.7), (2.13) over Goo, Xi, (p, V 7 ) V 7 gives the following system of equations 
(cf.(2.9)-(2.11) ; Goo =-1) 

N N 

c-]TA 2 +0 2 = i [iV'l' + E m?e- 2 ^ - 4) |^| 2 
i=i i=i 



+m +C£™y K -±m + 2 ^] , ( 2 - 14 ) 

i=l 

A- = |(m 2 e- 2 ^H 2 -m 2 e 2A *H 2 ) , (2.15) 

^-E a2 = i(h 2 + h ) > (2- 16 ) 



i=l 
N 



1? - ^ + (J2 rfe~ 2X * ~ 4)V + ^ = , (2.17) 



i=i 

N 



dU 

Eqs.(2.14)-(2.18) are invariant under the duality transformation which is the combination 
of the transformation (2.6) of "moduli" and of the interchanging of the momentum and 
winding modes, 

Xi — > — X{ , ip — > if , ip <-> %j) , rrii <-> rhi . (2-19) 

The system (2.14)-(2.18) is consistent in the sense that the time derivative of the "zero 
energy" constraint (2.14) is proportional to eqs.(2.15)-(2.18). In fact, eqs.(2.14)-(2.18) 
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have the following "mechanical" interpretation. Eqs.(2.15), (2.16) (combined with (2.14)) 
and (2.17), (2.18) follow from the action (2.7), (2.13) (with G 00 = -1) 

S = - j dt{K - U) , (2.20) 

where the kinetic and potential terms are 

N _ 2 



i=i 

N N 



U = e~v [-c+±(Y, m i e ~ 2Xi -*M 2 + l(£™y Xi -*M • ( 2 - 21 ) 

i=l i=l 

The action (2.19)-(2.21) describes a mechanical system with the M = 1 + N + 4 dimen- 
sional configuration space with the Minkowski signature metric qab = ^Vab > Vab = 
(— , +, +) and the dilaton ip playing the role of a time coordinate . This space is curved. 
For example, its scalar curvature is R = \{M — 1)(M — 2)e~ ip ; it vanishes if N = 1 ( i.e. 
D = 2) and the matter fields are absent. Eq. (2.14) is then recognised as the zero energy 
condition K + U =0 which should be added "by hands" . 

It is interesting to note that eqs.(2.15), (2.17), (2.18) generalise the renormalisation 
group equations for the couplings which correspond to perturbations of the torus vacuum 
by marginal operators. Let us consider for simplicity the case of one dimensional torus 
(N = 1). The action which generates the corresponding (3 - functions is given by [20] 

L = J dxdx [ \G-\ d -^) 2 + IG- 1 ^? ~ 2T 2 + |T 3 + ...] , (2.22) 



G = G~ l , G 



c 



2A 



f?=T=\%, f = G=-\^. (2.23) 
The rotational invariance constraint = implies 

oo oo 



m= — oo 



so that the action and RG equations take the form [13,20] 

oo oo 

l = \y, ( G ~ lm2 - 4 ) i^-i 2 + 3 E ( Grh2 - 4 ) i^-i 2 + w 3 ' > ( 2 - 25 ) 

m=l m=l 

V'm = = 4(^ -1 ^ 2 - 4)V>m + 0(iP 2 ) , (2.26) 

<L = = \(Gm 2 - 4)^ + 0(^ 2 ) , (2.27) 



OO OO 

! (//. I 

G=- 



' t = -j(E w2 i^i 2 - g2 E^ 2 i^i 2 )- ( 2 - 28 ) 

m=l m=l 



2 dC" 1 ! 

Here ^ m , Vm ? C = e2A are "running" couplings , and the derivatives are taken with 
respect to the RG parameter r. The correspondence between eqs.(2.26)-(2.28) and (2.14)- 
(2.18) is established as follows. If we drop the second derivative terms in (2.15), (2.17), 
(2.18), set N = 1 , ip=A and keep just two modes ( m and in) then eq. (2.15) becomes 
identical to (2.28) while eqs. (2.17), (2.18) reduce to (2.26), (2.27). As discussed in [18] 
this prescription corresponds to taking a semiclassical limit if the cosmological equations 
(2.14)-(2.18) are interpreted as RG equations in the presence of quantized 2d gravity. To 
satisfy (2.14) one should have c + (p 2 ^ 0. Given the value of c in (2.1) (i.e. c = 16 for 
D = 2, a' = 1) this implies that the RG parameter r should be identified with with the 
euclidean time r = it. 

We shall return to the analysis of eqs.(2.14)-(2.18) in Sect. 3. 

3. Let us now consider the case when the matter action is represented by the quantum 
vacuum energy of the string modes, i.e. by the partition function of the string compactified 
on a torus. In general, the partition function has the following duality invariant structure 



oc 



z(\, ^) = J2 e2(n_1)0 z m = E e(n ~ iv /-( A ) » ( 2 - 29 ) 

n=l n=l 

Z(\i ,<p) = Z{-\ ,<p) , f n {\ ,<p) = f n (-Xi ,<p) . (2.30) 
In the large radii limit 

JV 



Z n (\i -> oo) = d n V , ]/ = exp^Ai , (2.31) 



i.e. Z is proportional to the volume factor. Comparing with field theory and assuming 
that string ultraviolet cutoff corresponds to a proper time cutoff one may say that it is 
the momentum mode contribution that controls the large radius limit. In that sense the 
contributions of the winding modes (that make Z symmetric under (2.6) ) control the 
small radii limit 

N 

f n (\i -> -oo) =d n exp(-nJ2K) , h(\ -> -oo) = d{V~ x . (2.32) 

This interpretation is, in fact, ambiguous since from the cutoff field theory point of view 
both momentum and winding modes contribute similar terms ( e.g. with positive and 
negative powers of radii) to the partition function. The temperature dependence can be 
included by compactifying the euclidean time [21] (so that the inverse temperature (3 will 
couple to \J — Goo ) an d does not change the duality invariance property of Z . Assuming 
adiabaticity of evolution we may replace \i,<p and (3 by functions of time and take the 
matter action in the form (see also [22]) 

S m = J dW-Goo F(\,ip,f3^-G 00 ) , F = VT , Z = —(3F . (2.33) 

Let us consider, for example, the one-loop contribution when F does not depend on <p 
( this is a good approximation if the effective coupling e v is small). The energy-momentum 
tensor can be represented in the form 

V = ^=^L = diag(p , e !1 >, , .... e M »p N ) , (2.34) 
- 2SS - E E =F + I 3 d -L, (2.35) 



V 5G oo V dp 

^ = -vlx-=v ' P * = ~dK ■ (2 - 36) 

The conservation of the energy (following from the invariance under reparametrisations of 
time ) implies (E = dE/dt ) 

N N 

p + J2*i(P + Pi) = Q , E + J2*i P i = V ■ ( 2 - 37 ) 

i=l i=l 



Since F = F(\(t), f3(t)) eq.(2.37) is equivalent to the conservation of the entropy S = 

l3 2 dF/d(3, E + Eti >«Pi = S/P • 

The resulting cosmological equations (2.9) - (2.11) are [12] 

N 

c-J2^ + ^ = ^E , (2.38) 

i=i 

\ i -<p\ i = ±e lfi P i , (2.39) 

N 

0i-J2 X2 i = 1 ^ E ■ ( 2 - 4 °) 

i=l 

The time derivative of (2.38) vanishes on eqs.(2.39), (2.40) if eq.(2.37) is satisfied. 
Eqs.(2.37)-(2.40) are invariant under the duality transformation (2.6) if F is duality sym- 
metric (note that under the duality E — > E , Pi — > —Pi ). 

Solving the adiabaticity condition one can in principle express the temperature in 
terms of Aj . Then eqs. (2.38)-(2.40) will take the same form as in the zero temperature 
case with E and Pi replaced by £(A) and —dS/dXi. 

Re- written in terms of the original dilaton field <j> (see (2.5)) eqs.(2.38)-(2.40) take the 

form 

N N 

c-^Af + (2^-^A,) 2 = e 2 ^ , (2.41) 



i=i i=i 

N 

A i -(20-X)A < )A i = |e 2 *p i , (2.42) 
i=i 

N N 

20-^A i -^A 2 = ie 2 ^p . (2.43) 
i=i i=i 

It is easy to see that if c = a solution with = const may exist only if the energy- 
momentum tensor (2.34) is traceless, i.e. 

N 

J> = p . (2.44) 

i=i 

The latter relation ( which is not duality symmetric) is not in general satisfied by the free 
energy which actually appears in string theory. This point was already emphasized in [12] 
(cf. [3,4]). 
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A simple way to see that <p = const is not in general a solution of the cosmological 
equations if the one-loop vacuum energy is used as a source is to start with the action 
rewritten in terms of the "rescaled" metric (cf.(2.1),(2.33) ) 

G'^ = e-^ D -^G^, (2.45) 

S = - J d D xV^G e~ 2 ^ [c+ R+ A{dcj)) 2 } + J d D xV^G T = 

- J d D xV^G' '[Rl - -p^(d<P) 2 + ce 4 <^ D - 2 ) ] + Jd D xV^G 7 e 2D ^^ D - 2 ^. (2.46) 

The equation for which follows from (5.10) is not solved by <j> =const if c = and F^O. 

4. Before turning to the analysis of eqs. (2.9)-(2.11) with duality symmetric matter 
let us recall their solution in the absence of matter [15]. If S m = eq.(2.10) implies 

A; = hi e^ , ki = const . 

Substituting this into eq.(2.9) we finally get ( assuming D > 26 ) 

<f = ip G — In sinh 2bt , = A^o + Qi In tanh bt , (2.47) 
N 1 

J2lf = 1 . «V =-(£»- 26) , q i = k i b~ 1 e^ . 
i=i 

The solution for D < 26 is obtained by the substitution b — > ib ; in particular , for D = 2 

V? = (fio - In sin 26* , A = A ± In tg bt , a'b 2 = 4 . (2.48) 
The solution for the zero central charge ( D = 26 ) is given by 

JV 

cp = (po - In t , Xi = Xio + qi In t , ^q 2 = 1 . (2.49) 

i=i 

This is a Kasner-type solution with some radii infinitely expanding and others infinitely 
contracting. Such behaviour will be changed by the presence of matter. Note that the 
duality transformation (2.6) corresponds to qi — > —qi , i.e. it relates contracting solutions 
to expanding ones. 
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3. Solutions with classical matter 

1. In this section we shall study the solutions of the system of equations (2.14)-(2.18). 
We shall make a number of simplifying assumptions. In what follows we shall consider the 
isotropic case of all time dependent radii being equal A^ = A. This assumption is consistent 
with eq.(2.15) if the modes ifj and ifj are such that rtii = m , fhi = m , i = 1, N. The 
coresponding metric (2.4) 

N 

ds 2 = -dt 2 + e 2X ^J2 dx2 i > 
can be also interpreted as a spacially flat isotropic cosmological metric. Then 

(p = 2(f) -NX , V = e NX . (3.2) 

To be able to change the value of the central charge we shall assume that there may be a 
number K of "static" spacial dimensions, i.e. 

D = l + N + K , c = --(N + K -25) . 

3 

In general, the duality transformation (2.6) 

A = —A , = (p 

corresponds to inverting the direction of evolution of the scale factor 
Eqs. (2.14)-(2.18) reduce to 

c - NX 2 + <p 2 = \ [ + (Nm 2 e~ 2X - 4)|^| 2 

+ |^| 2 + (A^m 2 e 2A -4)|^| 2 + 2t/] , 
A - 0X = \ ( m 2 e- 2X \i;\ 2 - m 2 e 2X \^\ 2 ) , 

V3 - iVA 2 = i ( H 2 + |^| 2 ) , 
12 



(3.3) 



(3.4) 
(3.5) 
(3.6) 



2„-2A ■ 9U 



ip-<pip + (Nm 2 e~ 2A - A)ip + 



dip* 



~ ~ ~ BU 

ifj-¥ifj + (Nm 2 e 2X - 4)V> + = . (3.7) 

dip* 

We shall first ignore the dynamics of the matter fields ip and ip and consider them as 
"static" sources. Such an assumption may be justified if the the full interaction potential 
( including U as well as interactions of ip and ip with other modes which we ignored) has 
extrema (see, however, the discussion at the end of this section). Then eqs. (3.4)-(3.6) 
take the form 

-NX 2 + ip 2 = 2NW(X) , W = \\x e" 2A + \ji e 2A - C/2N , (3.8) 

X - ipX = e~ 2X - fx e 2X = -W (X) , (3.9) 

<p-NX 2 = 0, (3.10) 

where 

^i = \m 2 \^\ 2 , li=\m 2 \^ , C = c + 2\ij\ 2 + 2\^\ 2 -U(ij^) . (3.11) 

Note that the derivative of eq.(3.8) is still vanishing if (3.10) and (3.11) are satisfied. Let 
us discuss some general properties of eqs. (3.8)-(3.11) assuming that c , m and m are such 
that W is strictly positive. Then eqs. (3.8), (3.10) imply that cp ^ and is always growing. 
Hence if we start with ip < (i.e. the effective coupling e v decreasing with time ) then (p 
will increase towards positive values but will never reach zero. Eq.(3.9) then describes a 
motion of a particle in a potential W > (which grows at large positive and negative A) 
with a damping term proportional to —p [12]. Introducing the energy of the particle £ we 
find from (3.8)-(3.10) 

£ = (pX 2 < , ip 2 = 2N£ > , (3.12) 

£ = lX 2 + W(X) , 
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i.e. 8 will decrease with time. The particle coordinate X(t) will thus be oscillating between 
its maximal positive and negative values with the amplitude of oscillations decreasing to 
zero (i.e. the particle trajectory be reflecting from the walls of the potential moving down 
towards its bottom). If at the initial moment the scale factor a = e A is increasing (A > 0) 
the expansion will continue untill a maximal value of a is reached. After the turning point 
the contraction will start untill a minimal value of a is reached , etc. 

The asymptotic values of A and ip correspond to the minimum of the potential W 

A* = \\n (n/p) , (p* =(p -Qt , 

Q=[2NW(X,)] 1/2 = [2N(^) 1 / 2 -C] 1 / 2 , W(X*) = 0. (3.13) 

If m = rh , ifj = if; then A* = and duality is asymptotically restored at t — > oo. To 
find the behaviour of the solution at large t we can approximate W near A = A* by the 
oscillator potential 

W = A + \Bi 2 , A = A* + £(£), A = Q 2 /2N, B = (Q 2 + C)/N = 2(/jjl) 1/2 , 

£+Qi + B£ = 0, £~e^, u = -\(Q± [Q 2 - 4B] 1/2 ) . (3.14) 

Depending on the sign of Q 2 —AB = 2{N — 4)( / u/i) 1 / 2 — C we get either exponential falloff 
or oscillations with exponentially decreasing amplitude. 

The solution in the case (p > can be obtained by time reversal (see (3.8)-(3.10)). We 
conclude that e^ will be increasing and A will oscillate with the amplitude of oscillations 
exponentially growing with time. 

One can also analyse the solutions of (3.8)-(3.10) by studying the behaviour of the 
trajectories on the (A , <p) plane. We have 

/ = NX 2 , A 2 = 1( f 2 - 2NW ) , f = <p , (3.15) 

/' = Jl = ±y/(p-2NW)/N . (3.16) 
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(near the tuning point A = eq.(3.10) should be used as well ). The trajectories on the 
(A, /) plane will lie above and below the limiting curves / = ±y/2NW(X). Consider, for 
example the case when / > , A > , A>0. It is easy to see that the trajectory will 
move up ( / > ) towards the right branch of the limiting curve and will hit it in a finite 
time. We shall present the explicit solution of eq.(3.15) in the region of large A (and for 
N = 1 ) in Appendix A where we shall demonstrate that the growing X(t) always reaches 
a maximal value. Since at the turning point / = , A = but A < the trajectory will be 
reflected back : A will start decreasing , the trajectory will go up until it will finally hit the 
left branch of the limiting curve ( a will reach its minimal value). If we start with / < 
the trajectory will lie below the lower limiting curve. / will increase along the trajectory 
so the amplitude of oscillations of A will be decreasing asymptotically to zero. Again , the 
two cases ( / > and / < ) are related by time reversal. 

2. The above qualitative analysis of the system (3.8)-(3.10) can be repeated for the 
general system (3.4)-(3.7) with time dependent matter fields. Since according to (3.6) ip is 
again non-negative <p will remain negative if we start with decreasing ip. Then the damping 
terms in the equations for the matter fields (3.7) will cause ip , ip to decrease with time 
towards their "background" ( presumably constant) values. The "potential" W in (3.8) 
will now include the positive kinetic terms for ip and ip and hence the energy £ in (3.12) 
will still satisfy £ < 0. Since the derivatives of the matter fields will decrease with time 
the qualitative behaviour of the solution for A ( oscillations with decreasing amplitude ) 
will remain the same.0 

To summarise, if we start with the initial condition (p < , i.e. the effective coupling 
e v decreasing with time ( this is a natural assumption in order to justify the neglect of 
quantum corrections in the matter action ) then the scale a oscillates between the maximal 

1 If the structure of matter interactions is such that there is no non- vanishing backgrounds for 
ip and ip then the matter field contributions to eqs.(3.4) - (3.7) will be decreasing to zero. Hence 
the late time solution will be approaching the "free" solution (2.47) and the period of oscillations 
will become infinite. 
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and minimal values. If the matter action is "self- dual", i.e. \i = fx then a max and a m i n 
asymptotically approach the Planck scale \foi! = 1 , i.e. A = 0. If one starts with a large 
value of a it will take, of course, a large number of oscillations before a max will become of 
Planck order. 

As it is clear from the structure of the potential W the role of the winding mode 
contribution (~ exp2A) to the "classical" system of equations (3.4)-(3.6) is to prevent 
infinite expansion . The winding mode term in eq. (3.4) is thus opposing the expansion 
[12]. Without it we would have either the expansion to A = oo or first contraction to some 
minimal radius and then infinite expansion. The role of the momentum mode contribution 
is "dual": it prevents contraction to A = — oo , i.e. provides the existence of a non-zero 
minimal radius of contraction. 

One may question the validity of the assumption used in deriving (3.8)-(3.10) that 
the matter fields in (3.4)-(3.7) may be taken to be constant. It may be more appropriate 
to set constant the rescaled fields ft = exp(— if/2)ijj , ft = exp(—<f/2)ip for which there is 
no if term in the analog of eq.(3.7). Then eqs.(3.8)-(3.9) will have the same form but with 
extra factor of exp</? appearing in front of W (eq.(3.10) will be modified accordingly to 
preserve the consistency of the system (3.8)-(3.10) ; being redundant, it may be ignored). 
Similar system appears in the one-loop "quantum" case discussed in the next section and 
the behaviour of the solution is similar. Let us note that the presence of the exp ip factor 
in the r.h.s. of equations corresponding to the classical string modes is suggested also by 
the form of the combined action (2.7) where S m is given by the the classical string (a 
model) action. The absence of the dilaton factor in front of the classical string action 
implies its presence in the r.h.s. of the corresponding cosmological equations [10]. Hence a 
condensate of classical string modes is probably described by the rescaled classical fields. 
The presence of the dilaton factor in the r.h.s. of equations describing the cosmological 
evolution in the case when the "matter" is represented by classical string modes was also 
assumed in [12] where the conclusion that the winding modes oppose the expansion was 
already made. 
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4. Solutions with quantum matter 

1. In this section we shall study the cosmo logical evolution in the case when the 
matter action in (2.2)-(2.3) is represented by the (zero temperature) one loop partition 
function of the string compactified on the torus. The general form of the solution will 
be the same as in the case of the classical matter discussed in Sect. 3 : the scale factor 
will oscillate between its maximal and minimal values with the amplitude of oscillations 
decreasing to zero or increasing to infinity depending on whether the dilaton coupling is 
decreasing or increasing with time. 

As in Sect. 3 we shall consider the isotropic case, 

dF 1 dF 

A, = A , F(X 1 , A at) = F(X) , P % = -— = --— = P . (4.1) 

Then eqs.(2.38)-(2.40), (2.37) take the form 

c- NX 2 +ip 2 = e v E , (4.2) 

A-^A = ie^P , (4.3) 
(p - NX 2 = \^E , (4.4) 
E + NXP = . (4.5) 

2. For completeness let us first discuss the case of a non-zero temperature (F = 
F(X(t) : f3(t)) ) and consider (4.2)-(4.5) as a "phenomenological" system of equations which 
should be supplemented by some equation of state relating E and P ( for more general 
analysis of the finite temperature case see [12]). In the case of the usual radiation - type 
condition 

p = 7 p , P = 1 E , (4.6) 

we find from (4.5) 



E = E exp(-N-fX) , P = 7E0 exp(-A^ 7 A) 
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(4.7) 



Then a particular solution of (4.2)-(4.4) with c = is given by ( cf. (2.49) ; similar 
solutions were discussed in [22,11] ) 

<p = <po + s In t , A = Ao + q In t , (4-8) 
S = -2/(l + 7 2 iV) , Q = 2 7 /(l + 7 2 iV) , (4.9) 
E =4(l- 7 2 iV)(l + 7 2 A^)- 2 exp(iV7Ao - v?o) , s - -fNq = -2 . 
One should assume -f 2 N < 1 to have E > . If 7 = 1/N (so that (2.44) is satisfied) 

s = —2N/(1 + N) , ? = 2/(1 + AT), 0= I^ + jVA) = const . (4.10) 

This solution describes a power - law expansion with the dilaton coupling decreasing with 
time. 

The above equation of state (4.6) explicitly breaks the duality invariance of (4.2)- 
(4.5). As we shall see, the behaviour of the solution is very different in the case when the 
duality invariance is preserved. The simplest duality invariant equation of state is P = 
, i.e. E = Eq = const. This seems to be a natural choice for an equation of state near the 
Hagedorn temperature [2]. In general , if the temperature is approximately constant and 
F is duality symmetric such a relation will always be true near the Planck scale A = 0, i.e. 
E(X) = E(-X) ss const , P = -F'/N w . In that case we find from (4.2)-(4.4) 

A = Ae v , ip 2 = NA 2 e 2ip + E^ - c . (4.11) 

For c = [12] 



1 4E 1 t-2VNA/Eo 

V = hi — ^ , A = AnH == m = , (4.12 

Y l E 2 t 2 -4NA 21 VN U + 2VNA/E0 

i.e. we conclude that e A is asymtotically approaching its maximal or minimal values in the 
case of expansion ( A > ) or contraction (A < ) while the dilaton is decreasing with 
time. The duality transformation corresponds to A — > —A or, equivalently, t — > — t . Note 
that ip = const is not a solution of (4.2)-(4.4) (and 7^ for 7 = , cf.(4.10) ). 

18 



3. Let us now turn to the analysis of the solutions of eqs.(4.2)-(4.5) for the duality 
symmetric zero temperature case, i.e. 

E = F(X) = F(-X) , P = -F'/N , F(X -> ±00) -> d t e ±NX . (4.13) 

Let us first formally assume that E is positive (we shall consider the case of negative E 
later on) . Though the Casimir energy in field theory is usually negative, the bosonic string 
vacuum energy is in any case not well defined (divergent due to tachyon) but the aim of 
our analysis is only to illustrate some general properties of the system (4.2)-(4.4). 

Assuming (4.13) one finds that eq.(4.5) is satisfied automatically and (4.2)-(4.4) take 
the form 

— iVA 2 + ip 2 = 2NW , W(X, Lp) = [e ip F(X) - c ]/2N , (4.14) 

dW 

\-<p\ = -—, (4.15) 

Cp — NX 2 = NW — \c . (4.16) 

Since F(X) ( and hence W(X) ) has the minimum at A = there exists the "fixed point" 
self-dual solution 

A = , ip 2 = F e v - c , F = F(0) > , (4.17) 

i.e. if c > and (p < 

Lp = ip - 2 In I sin \yfc t\ , e^ = c/F 

Solutions with generic initial conditions will "spontaneously break" duality. 

A qualitative analysis of a general solution of (4.14)-(4.16) is similar to that of the 
"classical system" (3.8)-(3.10). Let us consider the case of <p> < 0. We note that since Cp 
is positive Lp will grow but will remain negative if it was negative at the initial moment, 
i.e. Lp will be monotonically decreasing. Treating <p(t) as a given function we can then 
interpret eq.(4.15) as the equation of motion for a particle with the damping term ~ — Lp 
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(4.18) 



and the time dependent potential W. The "energy" of the particle is decreasing with time 
(cf.(3.12)) 

£ = \\ 2 + W , £ = ip(X 2 + e ip F/2N) < , ip 2 = 2NS > . (4.19) 

If we draw the graph of the potential W(X, t) (which is exponentially rising at A — > ±oo and 
has a minimum at A = ) at several consequent moments of time it will be moving down 
towards the A axis. We shall assume that c is such that W > . Since £ is decreasing 
the trajectory of the particle will also be moving down reflecting from the walls of the 
potential. As a result, the radius e A will be oscillating from maximal to minimal values 
with the amplitude of oscillations asymptotically decreasing to zero, i.e. the solution will 
be approaching the self-dual point , X(t — > oo) = 0. The solution in the region near A = 
where F const was already analyzed above (see (4.11), (4.12)) : A is asymptotically 
approaching its constant value (zero) while e^ is decreasing to zero. 

If we start with ip > the direction of evolution is reversed: the amplitude of oscilla- 
tions of A grows with time, i.e. the maximal and minimal values of A are asymptotically 
approaching ±oo . 

It is possible to relate the analysis of the system (4.14)-(4.16) to that of the system 
(3.15)-(3.16) which appeared in the "classical" case. Introducing the new time coordinate 
dr = dte v / 2 and assuming c = we find from (4.14)-(4.16) 

NX' 2 = f - F{\) , 2/' = / 2 -F(A), / = ¥,' = ^, (4.20) 

| = iv / iV[/ 2 -F(A)]^. (4.21) 

Eq.(4.21) is similar to eq.(3.16) and so is the behaviour of the corresponding trajectories 
on the (A , /) plane. 

In Appendix B we shall present the explicit solution of eqs.(4.14)-(4.16) in the asymp- 
totic region of large A which demonstrates that a finite extremum value of A is always 
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reached in a finite time ( i.e. any trajectory on the (A, /) plane always hits a boundary 
curve ). 

Let us now consider the case when E in (4.13) is negative. Then it follows from (4.14) 
that A can never vanish (we are assuming that c is non-negative). As a result, there can 
be no turning points, i.e. the expansion goes to infinity and contraction goes to zero. 

To illustrate what happens when W is negative let us consider the special case of 
N = 1. If there are no extra static space dimensions, i.e. D = 2 , c = 16 then we may 
use the expression for the vacuum energy (2.29), (2.30) of a D = 2 string compactified on 
a circle which is known to all loop orders [23] 

f, = d 1 (e x + e" A ) , h = d 2 (e 2X + e~ 2X + 10/7) , ... , (4.22) 

where d\ is negative (see also [12]). Assuming again that the effective coupling e v is 
sufficiently small so that we can consider only the one-loop term F = f\ let us describe 
another way of solving (4.14)-(4.16) for A = 1 in the large radius limit. If A is large 
eqs.(4.14)-(4.16) can be represented in the form 

p - \p = , p(p-2X) = d t e p - c , p = ip + A , (4.23) 

i.e. 

p = Be x , -2p + p 2 = d x e p - c . (4.24) 
Integrating once the equation for p ( introducing y = exp( — |p) ) we find 

p 2 = q e p — dip e p — c , q = const . (4.25) 

As a result, if the radius ( which according to (4.24) is proportional to p ) is large and 
growing, it grows indefinitely. 
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5. Concluding remarks 

1. It useful to clarify the meaning of the duality invariance (2.19), (3.3) of our systems 
of equations (3.4)-(3.7) and (4.14)-(4.16). Let us consider, for example, (3.4)-(3.7) and 
assume for simplicity that m = fa. Then the duality invariance implies that for each 
solution (A , (p , ip , ip ) there exists another solution ( —A , ip , ip , ip ). If we consider 
the simplified system (3.8)-(3.10) and further assume that fj, = fx then the potential in 
(3.8) W = 2 / ucosh2A — C/2N will (like the potential in (4.14)) have the exact symmetry 
under A — > —A. Though there will exist the symmetric ("self-dual") solution (cf.(3.13), 
(4.17), (4.18)) A = , ip = <po — Qt in general duality will be "spontaneously broken" (by 
the initial condition A(0) ^ 0) on a time-dependent solution. As we have seen, however, 
the duality is asymptotically restored at late times (for a generic choice of initial conditions 
for A , A and <p < ) : the oscillating solution X(t) — > at t — > oo. For all the oscillating 
solutions the duality is true in the following "average" sense : the average value of A(£) 
is zero. Also, the behaviour of a given solution at large positive and large negative A is 
similar . This is in agreement with one's expectation that duality should be a symmetry 
between processes at large and small distances. 

2. Let us now discuss a model in which some of time dependent space dimensions 
are "uncompactified" , i.e. the vacuum energy depends on their scale factor in a trivial 
way. Let us consider the zero temperature case and set c = 0. Denoting the common scale 
of compact (noncompact) dimensions by A (A) and their number by (n) we find from 



eqs.(2.38)-(2.40) 



-nk 2 - 



NX 2 + v? 2 



(5.1) 



A — (pA 



(5.2) 



A — (p\ = 



IjV-V +nA F' , 



(5.3) 



(p - n'A 2 - NX 2 = ±e v+nA F 



(5.4) 
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Assuming F is positive a qualitative analysis of this system is similar to that of eqs.(4.14)- 
(4.16). Since (p > the dilaton continues to decrease monotonically , providing the damp- 
ing terms in the equations for A and A . Combining eqs.(5.1) , (5.2) and (5.3) one finds 
(cf. (B.1),(B.2)) 

A-<p = <p(A-<p) , A = tp + Ae ip . (5.5) 

As a result, one concludes that A first expands to its maximal value and then contracts 
to zero (note that there is no winding mode contribution which may protect A from con- 
traction to zero) while A oscillates between maximal and minimal values approaching the 
self-dual point A = 0. 

Since A in (5.1)-(5.4) corresponds to the non-compact dimensions it is more natural 
to introduce another "shifted" dilaton field as an effective coupling (cf. (2.5)) 

Tp = 2(p- NX , ip = Tp-nA . (5.6) 

To "absorbes" only the volume factor corresponding to compact dimensions and like ip 
is invariant under the duality A — > —A . The weak coupling regime corresponds to Tp 
decreasing with time. Tp in general will not be monotonic but will be decreasing at late 
times (if <p < 0). Let us consider the special solution with constant "internal" scale A = 0, 
i.e. A sitting at the minimum of F. Then (5.1), (5. 2) and (5.4) become identical to the 
system (B.l) (with N, A, di replaced by n,A, F(0) ) the solution of which is given by 
(B.8),(B.9),(B.12). The behaviour of Tp 

ip = ip + nA = ip + (y/n + l)/(Vn - 1) In z - 2(n + l)/(n - 1) ln(l + z) (5.7) 

is similar to that of A : it first grows with t, reaches its maximum (at z = (n + 1) 2 / (n — l) 2 ) 
and then starts decreasing. Since A reaches its maximum at a later time t (z = (n+l)/(n — 
1)) there exists the interval of t in which the non-compact dimensions are still expanding 
while the dilaton coupling is already decreasing. 
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3. Cosmological solutions for the space being a product of a number of large flat 
dimensions and an "internal" torus (and a finite temperature one-loop vacuum energy 
as a source) were previously discussed in refs.[3,4] with a conclusion that the presence of 
the winding modes prevents internal dimensions from expanding, stabilizing them at late 
times near the Planck scale. The dilaton dynamics was, however, ignored, i.e. the original 
dilaton <fi (see (2.5)) was implicitly assumed to be constant. We would like to emphasize 
that this assumption is inconsistent in general: as was already noted in Sect. 2 the system 
(5.1)-(5.4) does not have solutions with 

(j)= ^(if + nA + NX) = const 

unless (2.44) is true. The latter condition can be satisfied only in the limiting cases when 
the contributions of the winding modes can be ignored, i 

In this paper we considered only models at zero temperature. As we have already 
mentioned the more "realistic" finite temperature case was discussed in [12]. 

I would like to acknowledge J. Bagger and M. Tsypin for help and C. Vafa for col- 
laboration on ref.[12] on which most of the present paper is based. I am also grateful to 
Trinity College, Cambridge for a financial support. 

Appendix A . Asymptotic solution of the system (3.8)— (3.10) 

Here we shall present the solution of the "classical" system (3.8)-(3.10) in the case 
of N = 1. Our aim is to demonstrate the correctness of the qualitative analysis of the 
behaviour of the solution given in Sect. 3 which predicts the existence of a finite maximal 
radius of expansion. In order to solve eq.(3.16) explicitly we shall consider the region of 

2 Let us note also that no solution of (5.1)-(5.4) exists with Tp = const (the substitution of 
(p = _jv A into (5.1)-(5.4) gives F = ). 
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large positive A in which the potential W can be approximated by W = ^fte 2 . Assuming 
that A > and absorbing the positive constant ft into A we find that eq. (3.16) takes the 
form 

f = y r P^ , (A.l) 

or, equivalent ly, 

k' = -h(k), h = k - Vk 2 - 1 , fc(A) = /(A)e" A , k = \{h + h~ l ) . (A.2) 
Integrating over h we get 

-|ln \h\ - = A-A = A . (A3) 

Here the integration constant Ao should be large positive. Suppose that / > , i.e. 
/ > 1 , k > 1 so that h changes from 1 at k = 1 to at k = oo. Then A increases 
from some large negative value at initial large k to its maximal value —1/4 at k = 1. The 
maximal value of A ( or the value of Ao ) is determined by initial conditions for A and /. We 
have thus proved that for the given exponentially growing potential the rising trajectory 
on the (A, /) plane will always hit the boundary curve so that the maximal value of A will 
be reached in a finite time. 

The case of / < , A > is analysed in a similar way. Then h is strictly negative 
and (A. 3) is still valid. Both k and A increase from their initial negative values to their 
maximal values k = — 1 and A = —1/4. 

If A < at t = A will be decreasing towards negative values and will finally 
reach the region where the constant and e _2A terms in the potential W (3.8) cannot be 
ignored. The solution in asymptotic region of large negative A can be found by the duality 
transformation: A — > —A . The potential will be dominated by the term e~ 2A which will 
imply the existence of a minimal radius of contraction. The full solution can be found by 
"sewing" the solutions in the two asymptotic regions. 
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There is another way of representing the solution of (3.8)-(3.10) in the asymptotic 
regions of large ±A. Consider, for example, A < and A > , ip < , i.e. the case of 
expansion starting from small radius. Introducing the "conformal time" r 

dT = e~ x dt , F' = ^=e x F, (AA) 

(IT 

we find from (3.8)-(3.10) (N = 1) 

p" - p' 2 = , -p' 2 + 2pV = p , p = X + ip , (A.5) 

p = p - In |r| , (p = (p - iln \t\ - pr 2 /4 , r = r - r , (A.6) 
A = A - |ln |r| + pr 2 /4 , dt = drr" 1/2 exp(A + pr 2 /4) . (A.7) 

Appendix B. Asymptotic solution of the system (4.14)— (4.16) 

Below we shall solve eqs.(4.14)-(4.16) in the region where A is large and positive ( the 
solution for large negative A can then be found by the duality transformation). Then the 
potential W in (4.14 ) is given by (see (2. 31), (4. 13) ; we shall assume that d\ is positive) 

W = W/2N , W = d 1 e' p+NX , 

and our system takes the form ( the terms proportional to c can be neglected in the large 
A limit ) 

-NX 2 + cp 2 =W , \-(p\ = -±W , (p- NX 2 = \W . (BA) 
Combining these equations we get 

A-y3 = ^(A-<^), X = fi + Ae ip , (B.2) 

Cp = \{NX 2 +p> 2 ) , (B.3) 
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To have an expansion with decreasing dilaton , i.e. A > , ^<0we should take A < 
. Substituting the expression for A in (B.2) into (B.3) we get the equation for ip which is 
easy to solve by introducing the new time variable r 

dr = dte^ , (ff = ^ = e" > , X' = ip' + A , (5.4) 
dr 

so that 

ip" = \[N{ip' + Af-ip' 2 ] • (5.5) 
Let us first consider the case of N > 1 . Integrating (B.5) we find 

|4±^| = e -v^ , (5.6) 

-u = — = , v = —== . (5.7 ) 

As a result ( we assume that —u < <// < —v ) 

/ u + vz _ _ s /n At 

V? = — t~. , z = e 

1 + z 

¥> = ¥>o + - In * ~ ]\TTT m (1 + *) > ( S - 8 ) 

1 2 

A = A + = = In z-— — -lnfl + z) . (5.9) 

y/N(y/N-l) N — 1 K 1 



'N(VN - 1) 

The first equation in (B.2) 



cp' 2 -NX' 2 = d 1 e-^ +NX , (BIO) 



is then satisfied if 

A 2 = ^(iV-lH . (5.11) 
The relation between the original time t and z is 

dt = dre-^ = -(NA 2 )- 1 / 2 dzz a (l + z) b , a = ^2—, 6 = , (5.12) 



which can be approximated by t ~ z h = exp(\ZiV.A/iT) , < h < 1 , i.e. z decreasing to 
zero corresponds to t growing to infinity. 
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As it is clear from (B.9) , X(z) first grows to its maximal value at z* = and 
then decreases. The function (f(z) is monotonically decreasing with z — > . We conclude 
that a maximal radius of expansion is reached in a finite time ( z* corresponds to a turning 
point of the trajectory ). Once the contraction starts we eventually reach the region of 
small A where other terms in W (4.14) (which are present because of the contributions the 
winding modes) cannot be ignored . 

If N = 1 the integration of (B.5) gives 

(p = <po - \At + A- x e Ar , A = A + \At + A~ 1 e Ar . (5.13) 

A different way of solving (B.l) in the N = 1 case was discussed in Sect. 4 (see (4.23)- 
(4.25)). 
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